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Abstract 

Quantum state reconstruction involves measurement devices that are usually described by ide- 
alized models, but not known in full detail in experiments. For weak propagating microwaves, the 
detection process requires linear amplifiers which obscure the signal with random noise. Here, we 
introduce a theory which nevertheless allows one to use these devices for measuring all quadra- 
ture moments of propagating quantum microwaves based on cross-correlations from a dual-path 
amplification setup. Simultaneously, the detector noise properties are determined, allowing for 
tomography. We demonstrate the feasibility of our novel concept by proof-of-principle experiments 
with classical mixtures of weak coherent microwaves. 

PACS numbers: 03.65.Wj,07.57.Kp 
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The reconstruction of the Wigner function jl| or density matrix of a propagating quantum 
field represents a cornerstone in quantum optical measurement theory and experiments. In 
quantum homodyne tomography 2!, for example, the signal is combined with a local 
oscillator in a beam splitter, and the intensities at the output ports are subtracted to produce 
the measurement of the amplified field quadratures in terms of a histogram. The latter 
gives access to all quadrature moments, or, equivalently, the Wigner function In this 
procedure, it is of utmost importance that field amplifiers are not needed in the optical 
domain because photodetectors with sufficient efficiency are available at the single-photon 
level In contrast, in the 1 — 10 GHz range, which has become highly relevant due to 
the advent of circuit quantum electrodynamics (QED) jsUl, only theoretical proposals 

- Il4| . Consequently, 
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exist for the detection of propagating single microwave photons 
the detection of few-photon microwave signals requires linear amplification. Within well- 
established "off-the-shelf technology, cryogenic high electron mobility transistor (HEMT) 
amplifiers lend themselves to this purpose. They offer flat gain over a broad frequency range, 
but they obscure the signals by adding random noise 15|, ll6j of 10-20 photons at 5 GHz. 
Nevertheless, we prove here that it is still possible to measure all quadrature moments 
of few-photon propagating microwaves in this situation. Furthermore, we show that our 
proposed reconstruction method also produces a measurement of all quadrature moments 
of the detector noise. In this sense, moving from a single amplification chain to a dual-path 
configuration constitutes a step beyond pure state reconstruction and towards the complete 
calibration of the measurement device, i.e., detector tomography 17|. We note that so far 
only state reconstruction of the intra-cavity field has been demonstrated in circuit QED joj. 
However, quantum states of propagating microwaves themselves can be valuable in quantum 



information processing [18| and their full reconstruction represents an important open issue. 

We first develop a theory for the measurement of all moments of both a propagating 
quantum microwave signal and the noise added by the detector. The basic idea is illustrated 
in Fig. [IK. A signal S\/2 is equally split at low temperatures by means of a four-port 
50-50 microwave beam splitter. The outputs are amplified and synchronously digitized. 
During this process, the amplifiers add the independent noise contributions Xi ^"^^ Xi fo 
the s plit signals in the detection chains 1 and 2, respectively. Assuming a 180° hybrid 



ring 



19| as a beam splitter (cf. appendix), the recorded time traces are C\ = G{S^V\x\) 



and C2 = G'(— 5'+y+X2)5 where G is the gain and V represents the ancilla state present in 
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Figure 1: (a) Schematics of the dual path amphfication and detection setup. The ancilla port of 
the 180° hybrid ring is terminated by a 50$! load at T = 300 mK. The split signals are amplified 
and detected in separate chains, which add the noise xi and X2 to the signal, (b) Amplification and 
detection chain. The noise temperature is dominated by the cryogenic amplifiers anchored at 4.2 K. 
The isolation of the hybrid ring (40 dB) and the circulator prevents correlated noise contributions 
from the two chains. The phase difference between Ci and C2 is close to 180°. 



the fourth port of the hybrid ring. Since we assume full knowledge of V^, an obvious choice is 
the vacuum or, as in this work, a weak thermal state. Vacuum or thermal states at gigahertz 



frequencies can be prepared reliably in an experiment by controlling the temperature 2^, 21 |. 
In the next step, we compute suitable correlations of the form (CfC™), where £, mGNg 
and the brackets denote ensemble averaging. In contrast to Hanbury Brown and Twiss 
experiments based on intensity (power) correlations 20|, our method is devoted to the 
correlations of field quadratures (voltages). For the first signal moment, the mean values 
(Xi) and (X2) vanish and one obtains (S) = (Ci) /G = — (C2) / G. All higher moments of signal 
and noise can now be calculated by induction: 
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We note that our conjectures of equal gain in both chains and a perfect 180° hybrid ring are 
not a restriction in practice. In the derivation of the above formulas, the mutual statistical 
independence of S, V, Xi and X2 is crucial because it implies {S^V^XiXl) = i^^) (^'^) (Xi) (X2) 
for /3,7,5, eGNo. The latter formula also shows that Eqs. dl])-© are suitable for quantum 
signals, where S, V, xi, and X2 have to be interpreted as operators. 

Explicit expressions for moments of signal and detector noise can be calculated in the 
spirit of Eqs. (II])-(l3]). We assume {V'^^'^^) = ( j G N) for the ancilla as, e.g., in the case 
of Gaussian statistics. We note that for the initial correlation (CfC™) other choices than 
i = n—l and m = l are possible, as long as £+m = n and £, mGN. Typically, balanced 
products with ip^m result in simpler expressions because of higher symmetry. Starting from 



G{S ) = iC) = -{C2) 
G^{S^)= - {C^C2) + G^V) 

G'{S') = - {GIG2) - (Ci) {{Gl) + (C1C2) - 3G'{V')) 
G^{S^} = - G^{V^) - 6G^V^){GiG2) + QG^V^ 
+ {G,G2r + {G!G'2)-{G',){G'2) 

(Xi) = 

G'{xl) = {Gl) + {G^G2)-2G'{V') 
G'ixl) = (C?) + {C!C2) - 2{G,) {{GD + 



{GlCD we obtain: 
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G\x{) = {Ct) - l2G'{V^){CiC2) - 12G^{V^){Cl) 
+ Q{C^C2){Cl) + 5{C,C2f + 12G^{VY 
- MC^){C!) - 4{C,){GlG2) + 8{G^f{G',) 
+ 8{G,f{G,G2) - {ClGl) + {C\){Gl) 

Similar formulas can be derived for x^. Furthermore, we assume the first moment of the 
detector noise to vanish for both chains throughout the experimental part of this work. 
Practically, this implies an offset correction (c.f. appendix). Finally, from Eqs. ([I])-©, the 
central moments can be retrieved with the binomial transformation 



((5-(5)r) = E0' 



(4) 

fc=0 

We emphasize the practical relevance of the above theory by conducting proof-of-principle 
experiments with weak classical microwaves. The setup is shown in Fig. [1] As signals, we 
use pulsed coherent microwaves with a frequency of 5.85 GHz generated by a source at room 
temperature. A series of cold attenuators ensures that the thermal noise at the signal port of 
the hybrid ring is restricted to that of an effective 50 f2-termination at the base temperature 
of 300 mK. The source power at the input of the hybrid ring. Pin, is related to an equivalent 
number of signal photons on average (poa) as described in the appendix. Figure [2^ shows 
the ensemble average of a typical signal used in our experiments. The pulse duration of 



1 [xs mimics standard cavity decay times in circuit QED experiments |10l . |22| . We first 
demonstrate the suppression of the amplifier noise via cross-correlations. The auto- variance 
(Cf) — (Ci)^ of the ensemble is depicted in Fig. [2]d, where one immediately notices the large 
offset of 35.7xlO~^V^ due to the amplifier noise. In the cross-variance (C1C2) — (Ci)(C2), 
this offset is efficiently suppressed by two orders of magnitude. As expected for a coherent 
signal, the variances are flat and do not allow to distinguish between the "on"- and "off"- 
regions of the pulses. The fluctuations of the variance signals are smaller for the cross- 
correlation than for the auto-correlation by a factor of 1.6, see Fig. |2j:;. Next, we prove that 
our method works efficiently at the quantum level, i.e., for signals of few photons on average. 
To this end, we investigate the resolution limits of the constituents of the variance, mean 
value and cross-product. In Fig. [2li, the root mean square power inside the pulse region is 
plotted against the signal power at the input of the hybrid ring. We find a large dynamic 
range of the mean value extending over six decades down to 0.001 poa. This means that 
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Figure 2: Detection of coherent microwave probing pulses (10^ traces averaged): (a) Mean value of 
down converted signal (10 MHz), (b) Auto- (blue) and cross-variance (red), (c) Fluctuations of the 
auto- and cross-variance. The traces are obtained by subtracting the respective time averages, and 
the auto- variance trace has been shifted by 8x 10~^ V^, for clarity, (d) Dynamic range of the mean 
value. Prms is the root mean square (rms) power at the digitizer. The numbers above the data 
points indicate the number of photons on average, (e) Dynamic range of the cross-product (C1C2). 
The arrows indicate the values obtained from the time traces shown in (f). The rms amplitude 
inside the pulse region is calculated after subtracting the time average from each data point, (f) 
Power dependence of the cross-product time traces in the limit of few photons on average, (g) 
Discrete Fourier transform of the pulse region of the traces shown in (f). The arrows indicate the 
20 MHz peaks. 

pulse energies as low as 3.7x10"^^ J (0.01 poa) are still very well detectable. The power 
dependence of the cross-product is displayed in Figs. E^-g. Down to 2 poa, the pulse region 
is clearly visible (Fig. O). For Ipoa, a Fourier transform (Fig. |2|;) reveals that the signal 
component can still be identified. However, the associated peak has approximately the same 
amplitude as the largest noise peak in the spectrum. Hence, the detection limit of the 
cross-product (see also Fig. [2^) and therefore the one of the cross- variance is 1-2 poa. We 
note that this is not a fundamental limit, but is rather due to technical issues such as the 
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Figure 3: Statistical mixtures of phase-shifted pulses: (a) Mean value and cross- variance. The 
inset shows the pulse distribution for ip = 9Q° . The amplitude a is extracted with a numerical fit. 
(b) Cross- variance amplitude a plotted versus ip. Squares, 100 poa. Triangles, lOpoa. The solid 
lines are fits to the data, (c) Third central moment at 100 poa measured with a single amplification 
chain. Note that, as expected, the trace is not a simple sine function, but oc sin'^(a;t). The insets 
show the histograms: equally distributed (brown) and with finite skew (violet), (d) Third central 
moment measured with the dual path method. The dataset is the same as in (c). 

ensemble size, filter bandwidth, or bit resolution of the digitizer card. 

So far, we have studied time-independent ensemble variances because (/(i)^) — (/(^))^ = 
always holds for deterministic signals. Time-dependent variance signals require statistical 
mixtures of signals. As a first example, a sequence of coherent microwave pulses with 
alternating phase shifts ±(/9 is applied to the input line. This corresponds to a statistical 
mixture with an equally distributed histogram. The mean value and cross-variance are 
Acos{ip)sm{(jjt) and — acos^(ci;t), respectively, where a = A^siv?{^p) is the cross-variance 
amplitude and u;/27r = 10 MHz is the signal frequency after the mixers. The frequency 
doubling cos^(ajt)oc[l+ cos(2a;t)] is confirmed by the data shown in Fig. [3^1. Repeating the 
experiment for different phase shifts, we reproduce the expected sin^((y9)-dependence of a for 
10 and 100 poa (cf. Fig.Eb). 

Studying signal mixtures also allows us to address time- dependent third central moments 
{{S—{S)Y) and {{xi~ ixi))^) ■ They are proportional to the skewness of their associated 
probability histograms and vanish for Gaussian distributions. In particular, both quantities 
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are zero for an equally distributed histogram. Hence, to observe an oscillating third central 
moment a statistical mixture with an asymmetric histogram (cf. inset of Fig. [3t) is required. 
In the experiment, we again create a train of pulses with 100 poa, but this time 75% of 
them are shifted by the phase (f and 25% by —(p. Figure [3t shows the third central moment 
calculated from the data of a single amplification chain. In contrast to the case of an equally 
distributed mixture, a clear oscillating signal is visible within the pulse duration for a skewed 
histogram. However, there is an offset of approximately —7 xlO~^ V^. In Fig. [Sji, one can 
see that this offset disappears when also taking into account the data from the second chain 
and applying the dual-path detection scheme described by Eqs. ([I])-(11D. Following these 
equations, we can write down the third central moment of Ci as the sum of the third central 
moments of signal, noise, and ancilla state: 

((Ci - {C,)f) = G'iiS - {S))') + G'iixi - ixi))') 
+ G'{{V-{V))'). 

Since {(y—{V))^) =0 due to the the Gaussian statistics of V, it becomes obvious that the 
offset in Fig. [3b must be G'^ {{xi — {xi))"^) ■ In this sense, the noise of the detection chain shows 
non-Gaussian statistics. Although the exact origin of the latter inside the detection chain 
remains unclear, our experiments provide a first confirmation that the dual-path method is 
indeed capable of characterizing signal and detector noise moments simultaneously. 

In conclusion, we present a method based on off-the-shelf technology and signal recov- 
ery techniques to gain access to arbitrary moments of both weak propagating quantum 
microwaves and the detector noise. In particular, we allow for the use of linear amplifiers 
adding 10-20 noise photons to the signals. We find that it is crucial to move from a single 
to two amplification and detection chains and successfully perform proof-of-principle ex- 
periments with statistical mixtures of weak pulsed coherent microwaves. We demonstrate 
sufficient sensitivity for the first two moments and observe indications for a non-Gaussian 
statistics of the detector noise. The experiments indicate that our dual-path method is 
a suitable tool for detecting propagating quantum signals such as squeezed states from a 



Josephson parametric amplifier 



23 



261], Fock states 9|] leaking out of an on-chip resonator 



or non-classical microwave field states generated in a two-resonator circuit QED setup 



5| 
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Appendix 



The ancilla state is assumed to be well known and, hence, its quadrature moments can not 
be detected by our method. This is not a restriction in practice because vacuum or thermal 
states at gigal 



temperature 



lertz frequencies can be prepared reliably in an experiment by controlling the 
Furthermore, it is more relevant to characterize quantum states such 



as Fock states, where part of the Wigner function is negative. 

Beam splitter types for the dual-path method. Although in this work we always assume a 
180° hybrid ring as beam splitter, other choices are possible, as long as they provide enough 
isolation between the output ports. However, one has to keep in mind that all lossless. 



matched and reciprocal beam splitters must be four port devices 
the case of the Wilkinson power divider 



19 



19| . Sometimes, as in 



2l| . the fourth port may be hidden internally 



and therefore complicate the analytical treatment of the ancilla state. Nevertheless, the 
equivalent of Eqs. (l)-(3) of the supplemented paper can in principle be calculated for any 
beam sphtter. 

Offset correction. The raw data is divided into segments of 4128 traces equivalent to 0.5 s 
measurement time. For each of these segments, the time average is subtracted from each 
data point before any other manipulation. Effectively, this procedure acts as a high-pass 
filter eliminating slow drifts in the data. 

The calibration of the input line can be performed in situ in our setup. To this end, the 
50 Vt load at the ancilla port of the hybrid ring is temperature controlled using a heater and a 
thermometer. Except for small deviations, we can assume the hybrid ring, ancilla load, and 
effective 50 Vt load at the signal port to have the same temperature. Both loads inject thermal 
voltage fluctuations into the hybrid ring. The auto- variance of these fluctuations follows the 
well-known Planck function. Knowing the bandwidth of our setup from measurements with a 
spectrum analyzer (51 MHz), the gains of the ampliflcation and detection chains are inferred 
from a numerical flt of the Planck function to the auto- variance data 2]J]. We flnd gains 
of approximately 110 dB for both chains, where small residual gain asymmetries have been 
absorbed in a compensation factor already. Together with the 0.5 dB loss of the hybrid ring 
and the total transmission from source to digitizer, we can extract an input line attenuation 
of 94 dB. From this, we determine the signal power Pin at the input of the hybrid ring. 
The corresponding number of photons on average is the pulse energy PinT'puisc divided by 
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the energy quantum /ix5.85 GHz. Here, h is the Planck constant and Tpuisc = 1 (xs the pulse 
duration. 
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